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We unify two recent results concerning equilibration in quantum theory. We first generalise a proof 
of Reimann [PRL 101,190403 (2008)], that the expectation value of 'realistic' quantum observables 
will equilibrate under very general conditions, and discuss its implications for the equilibration of 
quantum systems. We then use this to re-derive an independent result of Linden et. al. [PRE 
79, 061103 (2009)], showing that small subsystems generically evolve to an approximately static 
equilibrium state. Finally, we consider subspaces in which all initial states effectively equilibrate to 
the same state. 
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INTRODUCTION 

Recently there has been significant progress in under- 
standing the foundations of statistical mechanics, based 
on fundamentally quantum arguments jl- 11|. In particu- 
lar, Reimann [HHI has shown that the expectation value 
of any 'realistic' quantum observable will equilibrate to 
an approximately static value, given very weak assump- 
tions about the Hamiltonian and initial state. Interest- 
ingly, the same assumptions were used independently by 
Linden et al. [H 0], to prove that any small quantum 
subsystem will evolve to an approximately static equilib- 
rium state (such that even 'unrealistic' observables on the 
subsystem equilibrate) . In this paper we unify these two 
results, by deriving the central result of Linden et aZ.Q 
from a generalisation of Reimann's result. We also offer 
a further discussion and extension of Reimann's results, 
showing that systems will appear to equilibrate with re- 
spect to all reasonable experimental capabilities. Finally, 
we identify subspaces of initial states which equilibrate 
to the same state. 



that all subsystems physically interact with each other. 
In particular, given any decomposition of the system into 
two subsystems H = Ha <£)Hb, equation {1} will not be 
satisfied by any Hamiltonian of the form H = Ha ®Ib + 
I a <8> Hb (unless either Ha or Hb is proportional to the 
identity) [H]. 

Theorem 1 (Generalisation of Reimann's result [l[). 
Consider a d- dimensional quantum system evolving un- 
der a Hamiltonian H = ^2 n E n P n , where P n is the pro- 
jector onto the eigenspace with energy E n . Denote the 
system's density operator by pit), and its time-averaged 
state by uj = (p(t)) t . If H has non-degenerate energy 
gaps, then for any operator A, 



* A ^(\tr(Ap(t))-tr(Au)\ 2 ) 
where \\A\\ is the standard operator norm fl2j . 
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EQUILIBRATION OF EXPECTATION VALUES. 

We prove below a generalisation of Reimann's result 
that the expectation value of any operator will almost 
always be close to that of the equilibrium state [l|. We 
extend his results to include non-Hermitian operators 
(which we will use later to prove equilibration of subsys- 
tems), correct a subtle mistake made in Q when consid- 
ering degenerate Hamiltonians, and improve the bound 
obtained by a factor of 4. As in @, 0|, we make one as- 
sumption in the proof, which is that the Hamiltonian has 
non- degenerate energy gaps. This means that given any 
four energy eigenvalues Ek , E\ , E m and E n , 

(Ek = Ei and E m = E n ) 
E k -Ek = E m -E n ^ ^ or (1) 

[E k = E m and E l =E n ). 

Note that this definition allows degenerate energy levels, 
which may arise due to symmetries. However, it ensures 
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This bound will be most significant when the number 
of different energies incorporated in the state, charac- 
terised by the effective dimension d c s , is very large. Note 
that 1 < d c g < d, and that rf c ff = N when a measurement 
of H would yield N different energies with equal proba- 
bility. For pure states d c ff = tr(w 2 ) -1 as in but it 
may be smaller for mixed states when the Hamiltonian 
is degenerate. 

The quantity A(^4) gives the range of eigenvalues when 
A is Hermitian, and gives a slightly tighter bound than 
the operator norm. Following (2j . we could improve 
the bound further by replacing A.(A) with a state- and 
Hamiltonian-dependent term [15j j . however we omit this 
step here for simplicity. 

Proof: To avoid some difficulties which arise when 
considering degenerate Hamiltonians, we initially con- 
sider a pure state p(t) = \yj(t))('ip(t)\, then extend the 
results to mixed states via purification. 
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We can always choose an energy eigenbasis such 
that \i(>(i)) has non-zero overlap with only a single en- 
ergy eigenstate \n) of each distinct energy, by includ- 
ing states \n) = P„|^(0)}/V(V>(0)|P n |V(0)} whenever 
(ip(0)\P n \ip(0)) > 0. The state at time t is then given 
by 



D -iE n t/h\ 



(5) 



non-degenerate energy gaps whenever H does. The ex- 
pectation value of any operator A for p(t) will be the 
same as the expectation value of A' — A eg) / on the total 
system, and we also obtain A(A') = A(A), \\A\\ = \\A'\\, 
and d' eS = d e fi. However, note that tr(a/ 2 ) does not equal 
tr(w 2 ). Using the result for pure states, we can obtain 
([2]) in the mixed state case from 



where c n = (n\i/>(0)}. This state will evolve in the sub- 
space spanned by {|n)} as if it were acted on by the 
non-degenerate Hamiltonian H' = E n \n)(n\. For any 
operator A, it follows that 

4= (|tr(A[p(t)-a,])| 2 ). 
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= ||A|| 2 tr (]>>„| 2 |n)(n| 
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In the fourth line, we have used the assumption that 
the Hamiltonian has non-degenerate energy gaps, in the 
sixth line we have used the Cauchy-Schwartz inequal- 
ity for operators with scalar product tv(A^B) and the 
cyclic symmetry of the trace, and in the seventh line we 
have used the fact that for positive operators P and Q, 
tr(PQ) < ||P|| tr(Q). This gives the weaker bound in the 
theorem. 

To obtain the tighter bound, we note that a a is in- 
variant if A is replaced by A = A — cl for any complex 
c. Performing this substitution with c chosen so as to 
minimize ||A|| we can replace ||A|| with ||A|| = A(A)/2. 

An extension to mixed states can be obtained via pu- 
rification, following the approach discussed in 0] . Given 
any initial state p(0) on H, we can always define a pure 
state |0(O)) on % <g> % such that the reduced state of 
the first system is p(0). By evolving \4>(t)) under the 
joint Hamiltonian H' = H ® / , we will recover the cor- 
rect evolution p{t) of the first system, and H' will have 
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This completes the proof. □ 

In Rcimann proves that a\ < A(yl) 2 tr(a; 2 ) when 
A is Hermitian and the Hamiltonian has non-degenerate 
levels as well as non-degenerate gaps. However, it ap- 
pears that there is a subtle mistake in Q when extend- 
ing this proof to degenerate Hamiltonians. Specifically, 
the step from equation (D.ll) to (D.12) in [2[ does not 
follow if the state has support on more than one energy 
eigenstate in a degenerate subspace. A counterexample is 
provided by the mixed state p(0) = r|0)(0|(g)/, of a qubit 
and a /c-dimensional system, with H — (|0)(1| + 11)(0|) <8>/ 
and A = (|0)(0| - |1)(1|) (g) I. In this case <j\ = ±, 
A(A) = 2 and tr(w 2 ) = ^, giving a\ > A(A) 2 tr(w 2 ) 
when k > 4. However, subsequently in [2j, tr(w 2 ) is re- 
placed by an upper bound of max„ tr(p(0)P n ), and this 
also upper bounds d~^, so later results are unaffected. 
Note that the bound given by Theorem 1 for the same 
example is satisfied tightly for all k, as d c g = 2 and thus 
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DISTINGUISHABILITY 

When A represents a physical observable and p(0) a 
realistic initial state, it is argued in [l[ that the difference 
between tx{Ap[t)) and tr(Aj) will almost always be less 
than realistic experimental precision. This is then taken 
to imply that p(t) will be indistinguishable from uj for 
the overwhelming majority of times. 

However, the fact that two states yield the same ex- 
pectation value for a measurement does not necessarily 
imply that they cannot be distinguished by it. For exam- 
ple, a measurement yielding an equal mixture of +1 and 
— 1 outcomes for one state and always yielding for a sec- 
ond state clearly can distinguish the two states, despite 
the expectation values in the two cases being identical. 
Furthermore, even though any particular realistic mea- 
surement cannot distinguish p(t) from to for almost all 
times, this does not imply that for almost all times, no 
realistic measurement can distinguish p(t) from oj. This 
is because the optimal measurement to distinguish the 
two states may change over time. Finally, the measure- 
ment precision is not easy to define for measurements 
with discrete outcomes. 

To address these issues, we first note that the most 
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general quantum measurement is not described by a Her- 
mitian operator, but by a positive operator valued mea- 
sure (POVM). For simplicity, we consider POVMs with 
a finite set of outcomes, which is reasonable for realis- 
tic measurements, as even continuous outputs such as 
pointer position cannot be determined or recorded with 
infinite precision [l6j . A general measurement M is de- 
scribed by giving a positive operator M r for each possible 
measurement result r, satisfying J2 r M r — !■ The prob- 
ability of obtaining result r when measuring M on p is 
given by tr (M r p). 

Suppose you are given an unknown quantum state, 
which is either p\ or p 2 with equal probability. Your 
maximum success probability in guessing which state you 
were given after performing the measurement M is 



(8) 



where 



D M ( Pl , P2 ) = |tr(M r pi) - tr(M r p 2 )|. (9) 



We refer to Dm(pi, P2) as the distinguishability of p\ 
and p 2 using the measurement M. Similarly, the distin- 
guishability of two states using any measurement from a 
set M is given by 



Note that 



Dm(pi,P2)= max D m (pi,P2)- 



< D M { Pl ,p 2 ) < D( Pl , P2 ) < 1, 



(10) 



(11) 



where D(pi,p 2 ) = ^ tr j pi — p 2 \ is the trace-distance, 
which is equal to Dm(pi, P2) when A4 includes all mea- 
surements. 



EFFECTIVE EQUILIBRATION OF LARGE 
SYSTEMS 



For typical macroscopic systems, the dimension of H 
will be incredibly large (e.g. For Avagardo's number Na 
of spin-i particles, we would have d > 10 10 " ), and it is 
unrealistic to be able to perform any measurement with 
this many outcomes, let alone all such measurements. For 
practical purposes, we are therefore restricted to some set 
of realistic physical measurements M.. In this case, we 
would expect M. to be a finite set, as all realistic experi- 
mental setups (including all settings of variable parame- 
ters) will be describable within a finite number of pages 
of text. 

We say that a state effectively equilibrates if 



This means that for almost all times, it is almost impossi- 
ble to distinguish the true state pit) from the equilibrium 
state to using any achievable measurement. 

We can obtain an upper bound on the average distin- 
guishability as a corollary of theorem 1. 

Corollary 1. Consider a quantum system evolving un- 
der a Hamiltonian with non- degenerate energy gaps. The 
average distinguishability of the system's state p{f) from 
oj, given a finite set of measurements M, satisfies 



(D M {p(f),u)) t < 



E M eME r A(M r ) < N(M) 



(13) 



where N(M) is the total number of outcomes for all mea- 
surements in M . 

The first bound will be tighter when measurements are 
imprecise, as each outcome is weighted by A(M r ) € [0, 1], 
reflecting its usefulness in distinguishing states (l7j . 

Proof: 



(D M ( P (t),cj)) t 



max D M{t )(p{t),u) 



< Y, (DM( P (t),u)) t 

= \ E E(i tr ( M ^))- tr ( M ^)i)< 

MEM r 



< 



< 



< 



MGM r 



N(M) 



□ 



(14) 



(D M (p(t),ij)) t « 1. 



(12) 



In realistic experiments, we would expect the bound 
on the right of (|13|) to be much smaller than 1, implying 
that the state of the system effectively equilibrates to ui. 
Consider again our system of Na spins. If d c s > d 01 , 
even if we take M. to include any experiment whose de- 
scription could be written in 10 19 words, each of which 
generates up to 10 21 bytes of data, we would still obtain 
(D M (p(t),Lo)) t <l/(10 w22 ). 



EQUILIBRATION OF SMALL SUBSYSTEMS 

Now consider that the system can be decomposed into 
two parts, a small subsystem of interest S, and the re- 
mainder of the system which we refer to as the bath B. 
Then H — Hs ®Hb, where Hs/b has dimension dg/B- 
It is helpful to define the reduced states of the subsystem 
p s (t) = tr B (p(t)) and u> s = tr s (w). 

In such cases, it was shown in that for sufficiently 
large d e g the subsystem's state fully equilibrates, such 
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that for almost all times, no measurement on the subsys- 
tem (even 'unrealistic' ones) can distinguish p(t) from co. 
In particular, when p(t) is pure and the Hamiltonian has 
non-degenerate energy levels as well as non-degenerate 
energy gaps, it is proven in 0] that 



(D(ps(t),ujs)} t < 




(15) 



Extending this result to degenerate Hamiltonians and ini- 
tially mixed states is discussed in 

We cannot recover this bound directly from (|13|) by 
considering the set of all measurements on the subsys- 
tem, because this set contains an infinite number of mea- 
surements. However, we can derive (|15p from Theorem 
1 by considering an orthonormal operator basis for the 
subsystem, given by the d 2 s operators 12 1 



F 



1 



(d s /c +fei) 



I 



~*s-\(l + k{)m.o&ds)(l\ (16) 



where ko,k\ <G {0,1,... ds — 1} and the states |Z) arc 
an arbitrary orthonormal basis for the subsystem. Then 
writing (ps(t) — us) = J2k ^k(t)F k we have 



(D(p S (t),LJ S )) t 




tr | | 



d s tr{Y,\ k (t)\f(t)F?F k ) 



kl 



kl 



(\ k (t)\Ut))MF?F k ) 




tr((p(t)-u>)F£®l)f 



(17) 



In the second line we have used a standard relation be- 
tween the 1- and 2-norm, and in the sixth line we have 
used Theorem 1 for the non-Hermitian operator F^® I. 
Note that \fd~s~Fk is unitary, and thus \\Fl <g) I\\ = -A=. 



UNIVERSALITY OF EQUILIBRIUM STATES 

We have so far been concerned with when states equi- 
librate, rather than the nature of their equilibrium state. 



However, one of the notable properties of equilibration is 
that many initial states effectively equilibrate to the same 
state, determined only by macroscopic properties such as 
temperature. Given a particular Hamiltonian and a set 
of realistic measurements M, we can construct a parti- 
tion of the Hilbert space into a direct sum of subspaces 
% = ® k %ki such that all states within Hk with large 
enough d c g effectively equilibrate to the same state fifc. 

One way to achieve this is to choose the subspaces 
such that each projector n^. onto Hk commutes with the 
Hamiltonian, and such that any two energy eigenstates 
in Hk are hard to distinguish, i.e. For some fixed e sat- 
isfying < e « 1, and all normalised energy eigenstates 

l<>,li>eft fc 



D M (\i)(i\,\j)(j\) < 



(18) 



When d e g is sufficiently large, it follows that all states in 
Hk effectively equilibrate to fifc = H k / tr(Uf.), as 

(D M ( P (t),n k )) t < (D M ( P (t),uj)) t + (D M (uj,n k )) t 

N(M) , x ^ (i\cj\i) 



tr(n fc )' 



(19) 



where the sums in the second line are over an eigenbasis 
of uj (which is also a basis of Hk), and we have used 
the fact that Dj^(p,a) satisfies the triangle inequality 
(D M (p, c) < D M (ft, t) + D M (t, a)) and convexity, 



Dm £ 



PiPi 



< 



2_,PiF>M(pi,v), 



(20) 



where pi > and pi = I. 

When Hk can be chosen to be a small band of energies, 
the equilibrium state Qk will be the usual microcanonical 
state. 



CONCLUSIONS 

To summarise, we have shown that two key results of 
0, |]| about the equilibration of large systems can be 
derived from very weak assumptions (non-degenerate en- 
ergy gaps, and sufficiently large <i e ff), and a single the- 
orem (Theorem 1). In particular, for almost all times, 
the state of an isolated quantum system will be indistin- 
guishable from its equilibrium state oj using any realistic 
experiment, and the state of a small subsystem will be 
indistinguishable from los using any experiment. 

Although the first result has a similar flavour to the 
classical equilibration of course-grained observables such 
as density and pressure, it is really much stronger, as it 
encompasses any measurement you could describe and 
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record the data from in a reasonable length of text, in- 
cluding microscopic measurements. The second result 
has no classical analogue, as it yields an essentially static 
description of the true micro-state of a subsystem, rather 
than the rapidly fluctuating dynamical equilibrium of 
particles in classical statistical mechanics. Given the dif- 
ficulty of proving similar results in the classical case, it 
seems that quantum theory offers a firmer foundation for 
statistical mechanics. 
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